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Abstract. We consider the scattering of gravitational waves by the weak gravitational fields of lens objects. We obtain the 
scattered gravitational waveform by treating the gravitational potential of the lens to first order, i.e. using the Bom approxima¬ 
tion. We find that the effect of scattering on the waveform is roughly given by the Schwarzschild radius of the lens divided by 
the wavelength of gravitational wave for a compact lens object. If the lenses are smoothly distributed, the effect of scattering 
is of the order of the convergence field k along the line of sight to the source. In the short wavelength limit, the amplitude is 
magnified by 1 -l- /<■, which is consistent with the result in weak gravitational lensing. 
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1. Introduction 

Ground-based laser interferometric detectors of gravitational 
waves such as LIGO, VIRGO, TAMA and GEO are currently 
in operation to search for astrophysical sources such as neutron 
star binaries, black hole binaries and supernovae (e.g. Cutler 
& Thorne E002t . The gravitational wave signals from these bi¬ 
naries are extracted from the data using matched filtering with 
a gravitational waveform template. If the gravitational waves 
pass near massive compact objects or pass through intervening 
inhomogeneous mass distribution, the gravitational waveform 
is changed due to the scattering (or the gravitational lensing) 
by the gravitational potential of these objects. The gravitational 
waves do not directly interact with matter (e.g. Thorne n'987> . 
but the gravitational lensing occurs in the same way as it does 
for electromagnetic waves. In this letter, we investigate the ef¬ 
fects of scattering by lens objects on the gravitational wave¬ 
form. 

In the gravitational lensing of light, the scattering is dis¬ 
cussed in terms of gravitational lensing under the geometrical 
optics approximation, which is valid because the wavelength 
is much smaller than the typical size of lens objects. But in 
the case of gravitational waves, since the wavelength is much 
larger than that of light, geometrical optics is not valid in some 
cases. If the wavelength is larger than the Schwarzschild radius 
of the lens, wave optics should be used (Peters 1X974! Ohanian 
119741 Bontz & Haugan [l9811 Thorne n'9831 Deguchi & Watson 
I1986> . This condition is rewritten as M ^ 10^Mo(//10^Hz) ', 
where lO^Hz is the typical frequency of gravitational waves 


for ground-based detectors. Hence, we use wave optics in this 
letter. 

In the geometrical optics for the lensing of light, the strong 
and weak lensing are distinguished by the convergence field 
K which is the ratio of surface density of lens S to a critical 
density Scr ~ c^lGD]^, where Dl is the distance to the lens 
(Kaiser fT9921 Bartelmann & Schneider I2001II . In the strong 
lensing regime, a- ^ 1, the multiple images of distant source are 
formed. But the strong lensing probability is small, ~ 0.1%, 
for high redshift sources. Hence, the weak lensing approxima¬ 
tion, K <si 1, is valid for most sources. We use the weak field 
approximation in the wave optics. 

In the past, Peters J 197411 studied the scattering by a point 
mass lens and a thin sheet of matter in the weak field approx¬ 
imation, and obtained the scattered waveform for these lens 
models. The gravitationally lensed waveform ( which is the so¬ 
lution of wave equation o ) was given in Schneider, Ehlers & 
Ealco ([1992), Sec.4.7 and 7, using the diffraction integral un¬ 
der the thin lens approximation. Recently, several authors have 
been studying wave optics in the gravitational lensing of grav¬ 
itational waves using this integral (Nakamura ll9981 Nakamura 
& Deguchi 119991 Ruffa flX^ Takahashi & Nakamura 120031 
Yamamoto l2003l Macauart l20041 Takahashi l2004l '). In this let¬ 
ter, we present another method to derive the solution of Eq.(Q 
in the weak-field limit. We treat the gravitational field of lens 
to first order, i.e. using the Born approximation, and discuss its 
validity. We use units of c = G = 1. 
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Fig. 1. Gravitational lens geometry for the source, the lens and 
the observer. The lens is distributed around the origin. The 
source position is rs = ( 0 , -Z)ls), while the observer posi¬ 
tion is To = (So,£>l) where So is a two-dimensional vector with 
|So| « £>l- and Dls are the distances from the lens to the ob¬ 
server and to the source, respectively. 0° is the incident wave, 
while (jp + (f)^ is the scattered wave. 

2. Gravitational Waves Propagating through Weak 
Gravitational Fields 

We consider the scattering of gravitational waves by the weak 
gravitational fields of lens objects. We consider the scattering 
of scalar waves, instead of gravitational waves, since the ba¬ 
sic equation is the same as the scalar field wave equation (see 
Sect. lie of Peters I1974t . The gravitational fields of the lens 
objects are described by the metric as ds^ - g^ydxf^dx" = 
- (1 H- 2U) dt^ + (I - 2U) dr^, where U(r) is the gravitational 
potential of the lens. Since the propagation equation of the 
scalar wave is ( ■\f^g^''dv(l>^ = 0, we have the wave equa¬ 
tion as 

( H- 0(r) = 4a/ t/(r)0(r), (1) 

where o) is the frequency of the gravitational waves and we set 
0(f,r) = 0(r)e‘"^ 

We show the lens geometry of the source, the lens and the 
observer in FigH The lens is distributed around the origin of 
the coordinate axes. The source position is rs = (0, -Z)ls), 
while the observer position is Tq = (So, £>l) where So is a two- 
dimensional vector with |So| ^ £>l- and Dls are the dis¬ 
tances from the lens to the observer and to the source, respec¬ 
tively. 

In the unlensed case t/ = 0 in Eq. Q, we write 0*’(r) as the 
solution of this equation. Including the effect of U to first order, 
the scattered wave at the observer is written as 

0(ro) = 0°(ro)+ 0^(ro), (2) 

where 0' represents the effect of scattering. Using Green’s 
function for the Helmholtz equation, 7 |r - r'|, we obtain 
from Eqs. Q and 0 as 

. .2 ^ Aco\ro-r'\ 

0*(ro) =- mV- --t/(r')0V)- (3) 

^ J |ro-r'| 


Note that the above result m can be used if the lenses are 
broadly distributed between the source and the observer since 
the thin lens approximation is not assumed. 

We use the spherical wave emitted by the source as (f/, then 
we have - Ae‘"''''''®7 |r-rs|. InEigHI 4 / represents the 
incident wave emitted by the source, while 4 / + 4 / represents 
the scattered wave. 


2.1. Geometrically thin lens 


We assume that the lens objects are locally distributed at the 
origin. Then we have |r'| <si Dl.ls and |So| ^ £>l in Eq. 0. 
Using the second-order Taylor series' for these small quantities 
r' and So, 4>^ is reduced to 


Ds 


J (4) 


0‘(ro) 

0O(ro) 7T DlDls 
where Ds - D]^ + Dls and we set r' = (s',z')- 

Here, is the geometrical time delay which is given by 
fd(s'. So) = Ds/( 2 Dl£»ls) X [s' - (Dls/T>s) So]^. Using the two- 
dimensional gravitational potential 4f(s') - 2 J dz' U(s', z'), the 
result in Eq.0 is reduced to 


^7ro) 


a> 


Ds 




(5) 


^“(ro) 27r Z)l£»ls 

The above equation is also derived by expanding if/ to first order 
in the diffraction integral (see Schneider et a/. ll992t . The sur¬ 
face density of the lens is defined as S(s) = f dzp(s, z), where 
p is the mass density. Since the potential 4j is written using the 
surface density as i/f(s') = 4 J^f^s"^s") In |s' - s"|, the result 
0 is rewritten as (Appendix A), 


^7ro) 

0°(ro) 




= 2iw V/S(s') 


Ei (iwfd(s',So)) - In 


T>ls 

-s 

Ds 


(6) 


where Ei is the exponential integral function : £’i(ijii:) = 
- dt e''lt (e.g. Abramowitz and Stegun n'970> . If we shift 
the potential i/r by a constant value 4 ^ 0 , in Eq. 0 is shifted 
by Hence we can choose i/^o so that the second term of 

Eq.0 vanishes, and we eliminate this term. 

Let us discuss the validity of the Born approximation. If 
the scattered wave 4 / + is not too different from the incident 
wave 4 /, this approximation is valid. We discuss the condition 
of ^ 1 for the two lens models, point mass lens and 

smoothly distributed lens, in the following subsections. 


2.1.1. Point mass lens 


The surface density is S(s) = M6^(s) where M is the lens mass. 
Then, in Eq.0 is rewritten as (see also Peters IT^7?I Sect. 


Ill A) 


^°(rn) 


= 2iMoj Ei 


I IcuDls 2 \ 

imps 


(7) 


' |r - r'l r - (r • r') jr + - (r • r')^ /Uj /2 for r » r'. 

^ In other words, this corresponds to an additional phase shift in the 
incident wave : 4>^ IjP (I - itJi^o). But, is not changed 

under this shift since 4>^4'o ~ O(U^) is negligible. Hence, we are free 
to choose 4/0 ■ 
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Fig. 2. Ei(\x) as a function of x. The solid line is the absolute 
value of El, the dashed line is the real part, and the dotted line 
is the imaginary part. 


In Fig|2 we show £’i(ijc) as a function of x. The solid line is 
the absolute value of E{, the dashed line is the real part, and the 
dotted line is the imaginary part. From this figure, except for 
small ;c 1, IFil is smaller than 1. Thus the Born approximation 
is valid 1 in Eq. 0 if the lens mass M is smaller than 

the wave length of the gravitational waves A - 2nla) except for 
small impact parameter Sq ^ 

Especially for large impact parameter Sg, since Ei{\x) ~ 
IX for X » 1, Eq. 0 is reduced to 


^“(ro) 



itU^LS 2\ 
2 DlZ)s 'V ’ 


( 8 ) 


where is is the Einstein radius on the observer plane : ie = 
(4 MDl-Ds/Z)ls)^''^- Hence if the observer position So is larger 
than the Einstein radius sg, the Born approximation is valid 
irrespective of the wavelength of the gravitational waves. 


2.1.2. Smoothly distributed lens 


We assume that the lenses are smoothly distributed on the 
z - 0 plane in Eig^ In the unlensed case, the wave propa¬ 
gates through s' = (Dls/T>s)So in the lens plane. Expanding 
the lens potential around it, 0' in Eq. 0 is rewritten as 


0Hro) 

0“(ro) 


—itOlf/ H- /<■-!- 


i DlDls 
4io Ds 


VIk + 0 


/ T>lT>ls \" 
\ wDs / 



.(9) 


The first term, being -iwi/'(s') at s' = (Dls/T>s)So, can be elim¬ 
inated for the same reason as the second term in Eq.0. In the 
second term, k is the convergence field along the line of sight 
to the source defined as a- = S(s')/Scr at s' = (Dls/T>s)So with 
2cr - Dsl (47 tDi^Di^s). The third term is a correction term aris¬ 
ing from the effect of finite wavelength. If the scale of the den¬ 
sity fluctuation ~ |VsA'/a:| is much larger than the Eresnel 
scale ~ or if the geometrical optics is valid 

(a; ^ oo), the third term is negligible. Then, the incident wave 
is magnified by 1 -H /c, which is consistent with the result in 


the weak gravitational lensing (e.g. Bartelmann & Schneider 
IWfl .If ^ 1, the Born approximation is valid. 


2.2. Geometrically thick lens 


We consider the lenses are broadly distributed between the 
source and the observer. It is easy to apply the previous result 
in Sect. 2.1 to this case. 0' in Eq. 0 is rewritten as 


0'(ro) 


Cclz’- 

fo) ^ J Jo • 


Ds 


(Ds - z') z! 
xU(s',z - Di^s) 


( 10 ) 


with fd(s', z'. So) - Ds/ [2z' (Ds - z')] x [s' - (z'/Ds)So]^. If the 
lenses have a thickness of Az in the z direction, the correction 
for this thickness in the scattered wave is of the order of 
Az/Ds from Eqs. 0 and (llOt . Thus, if Az <K Ds the thin lens 
approximation is valid. 

Especially, for smoothly distributed lenses the result 0 is 
valid but k is replaced by 


K — 4n 



dz' 


z'jPs-z') 

Ds 


pi(z /Ds)So,z' - Di^s), 


( 11 ) 


where p is the mass density of the lens. 


2.3. Two-point correlation function 


Recently, Macquart (120041 (hereafter M04) derived the correla¬ 
tion function in the wave amplitude of the two detectors under 
the thin lens approximation. He suggested that measurement 
of the correlation function provides the power spectrum of the 
mass density fluctuation. In this section, we derive it without 
the thin lens approximation, but within the limit of weak gravi¬ 
tational fields.^ We consider two observers at Tq andro-i-ri with 
|ril ^ Ifol- The mass fluctuation is usually characterized by the 
power spectrum defined as P(k) - J(p(r)p(r + r')}e'^''''d^r'. 
The correlation in the potential {/(r) and t/(r') is written as 

(U(r)U(r')} P(k) (12) 


We obtain the correlation in the wave amplitudes 0*(ro) and 
0'(ro + fi) from Eqs. 0 and (HU as 


<0Xro)^^*(ro + ri)> 




^iwlro-r'l 




To - r' 


^-iw|r„+ri-r"| 


To -H ri - r" 


,0*(r") g-ik'(r'-r") 


(13) 


It is useful to change the integral variables to x = r' - r" and 
y = (r' -H r") /2, and we assume that x is much smaller than the 
distances Di is. Then, expanding these small quantities x and 
ri in the exponentials, Eq. (O is rewritten as 

<0^(ro)^‘*(ro-Hri)) ^ J d^k-j^ P(k) j d^x j d^y 

X _i_ g-iw[(r~y+r'py) x+(r~y)'ri] ( 14 ) 

_|ro - yf Ts - yl^ 

^ The correlation function due to electromagnetic scattering was ex¬ 
actly obtained, not only for weak fluctuation but also for strong fluc¬ 
tuation. See references, Ishimaru ll978l Tatarskii & Zavorotnvi l 19801 
for a detailed discussion. 
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where "v denotes the unit vector for v. A is the amplitude of 
the incident wave (see sentences after Eq.Q)- Performing the 
integral in Eq. (US, we obtain 

<0'(ro)^'*(ro + ri)) = 16w^ |^°(ro)p JP(q) 

X Pr/z'e-'""--', (15) 

Jo 

where n = (-z'q/wDs, 1) and q is the (x,y) component of the 
wave vector k. In the geometrically thin lens at z = 0 plane, 
P{q) is replaced by P^^(q)6(z' - £> ls ) where P^° is 2D power 
spectrum. 

Let us consider a single power law model for the power 
spectrum as an example. This model can be used for the fluc¬ 
tuation of cold dark matter and gas (M04, Sec.5.2 (& 5.3). The 
power spectrum is P{k) - Pq (k/ko)^" for kmin < k < kmax, and 
P{k) = 0 otherwise. The index is n » 3 - 4, and k^nin kmax- 
Then, the integral in Ea. (ll5t is dominated by the fluctuation at 
the largest scale of l/kmin- The exact solution of integral J15t 
was given by the hypergeometric function, but we present an 
approximate solution for simplicity. Assuming that the separa¬ 
tion of two detectors |ri| is much smaller than the largest scale 
of fluctuation 1 /k^i„, we have 


<0'(ro)^'*(ro + ri)> = 16oj^ |^°(ro)f P(k™i„)k^„Dse-‘“ 


2 - n 6n 


, (16) 


where ri = (si,zi). The first term represents the dispersion of 
the scattered wave, while the second term is two-point correla¬ 
tion function."^ If the lenses have a finite thickness AL along the 
z-axis, Ds in Eq. (HI should be replaced by AL. 


3. Conclusions 

We have discussed the scattering of gravitational waves by the 
weak gravitational fields of lenses by using the Born approxi¬ 
mation. We consider the two lens models, the point mass lens 
and the smoothly distributed lens, and discuss the validity of 
the Born approximation. For the point mass lens, the effect of 
scattering is roughly given by the Schwarzschild radius M of 
lens divided by the wavelength 4. If M < T or if the impact 
parameter is larger than the Einstein radius, the approximation 
is valid. For the smoothly distributed lens, the effect of scat¬ 
tering is of the order of the convergence /c, and if /c « 1 the 
approximation is valid. In the short wavelength limit, the re¬ 
sult is consistent with the weak gravitational lensing. We de¬ 
rive the correction term due to the effect of finite wavelength 
in the magnification. The two point correlation function is also 
discussed following the recent paper (M04). 

Acknowledgements. We would like to thank the referee for useful 
comments to improve the manuscript. 

We comment on the first term. In the previous work (M04), he 
includes the correlation between the scattered wave being second or¬ 
der of potential OiifP) and the incident wave. Then, the correlation at 
= 0 is subtracted (see Eq.(23) in his paper), and in Eq. <14> 

is replaced by . Hence, the first term would vanish if we 

included the second order of ij/ in the scattered wave. 


Appendix A: Derivation of Eq.© 

Inserting if/is') - 4 j£/^s"2(s") In |s'- s"| into Eq.(|3, we 
change the integral variable from s' to u s s' - s". Then, we 
have 


0°(ro) 


DlDls 

X 


/ 




f 


duu\nu e‘“" Jo(J3u), 


(A.l) 


where a = wDs/(2 DlDls), = 2ci'|s" - (Dls/Ds)So|, and 
Jo is the 0-th order Bessel function. By using series JQijJu) — 
(n!)~^, the integral in Ea. dA.H is rewritten as. 


f 


du u\nu e‘“" Joifiu) - -[y H- In(-ia)] 

4a 


4a 4—i n! \ 4a/ k 

n=l ' ' Lk=l 


(A.2) 


where y is the Euler constant (e.g. Gradshteyn & Rvzhik l2000t . 
Using identity 2k=i ^l^-^Qdt[l-{\- f)"]/L we have 

z s (-£)■ z EXtC 


^ _e->/3V(4a) 


(A.3) 


In the second equality, we use £’i(ix) = y H- In (-ix) -i- 
2^1 (i-*)" / (n ■ n!)- Inserting Eas. dA.2> and (IA.3II into 
Ea. dA.H . we obtain Eq. 
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